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ABSTRACT 

We present new Spitzer IRS spectroscopy of Cygnus A, one of the most luminous radio sources 
in the local universe. Data on the inner 20" are combined with new reductions of MIPS and IRAC 
photometry as well as data from the literature to form a radio through mid-infrared spectral energy 
distribution (SED). This SED is then modeled as a combination of torus reprocessed active galactic 
nucleus (AGN) radiation, dust enshrouded starburst, and a synchrotron jet. This combination of 
physically motivated components successfully reproduces the observed emission over almost 5 dex in 
frequency. The bolometric AGN luminosity is found to be 10 12 L© (90% of Lm), with a clumpy 
AGN-heated dust medium extending to ~ 130 pc from the supermassive black hole. Evidence is 
seen for a break or cutoff in the core synchrotron emission. The associated population of relativistic 
electrons could in principle be responsible for some of the observed X-ray emission though the syn- 
chrotron self-Compton mechanism. The SED requires a cool dust component, consistent with dust 
reprocessed radiation from ongoing star formation. Star formation contributes at least 6 x 10 10 Lq 
to the bolometric output of Cygnus A, corresponding to a star formation rate of ~ 10 Mq yr _1 . 
Subject headings: galaxies: active - galaxies: individual (Cygnus A) - galaxies: jets 



1. INTRODUCTION 

As one of the nearest powerful radio-loud active galac- 
tic nuclei (AGN), Cygnus A provides an excellent lab- 
oratory to study the environment and activity of pow- 
erful AGN. The luminosity of the AGN in Cygnus A 
(which is predominantly expressed in th e infrared) is high 
enou gh to classify it as quasar (e.g., iDiorgovski et al.l 
[1991 . While not seen in total intensity , broad Ha is 
seen in polarized light ()Ogle et al.l 119971 ). lending sup- 
port for the existence of an obscured broad-line region 
(BLR). The polarized broa d lines were d e tected within 
the ionization cone seen bv I Jackson et al.l (1996), giving 
them a possible scattering origin in the NLR. 

The host galaxy of Cygnus A is a cD elliptical near 
the center of a cluster which appears to be undergo- 
ing a merger with clust e r of s imilar size (|Ledlow et al.l 
120051) . iTadhunter et al.l (|2003[) obtained Hubble Space 
Telescope (HST) Space Telescope Imaging Spectrograph 
(STIS) spectra of the nuclear region. Stepping the slit 
across the nucleus, a velocity gradient indicative of ro- 
tation around the radio axis was observed. Model- 
ing the velocity as due to the potential of a super- 
massive black hole (SMBH) and stellar mass distribu- 
tion (measured from a 1.6/zm NICMOS image) gives a 
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SMBH mass of 2.5 ± 0.7 x 10 9 M©. This measurement 
is con sistent with black-hole - mass-host-galaxy relations 
(e.g.. iMagorrian et all 119981: iFerrarese fe Merrittl 120001: 
IGebhardt et al.H2000t IGultekin et al. I12OO90- 

The large scale radio morphology shows prominent 
hotspots, lobes, and a radio core. Both a jet and 
counterjet are visible in Very Large Array (VLA) and 
Very Long Baseline Array (VLB A) observations (e.g., 
iSorathia et al.llT996l ). Due to the edge-brightened mor- 
phology it is classified as an FR II source (Fanaroff- Riley 
Type II; iFanaroff fe Rlie^lr97l . 

X-ray observations of Cygnus A are consistent with 
the presence of a hidden quas ar. Chandra ACIS obser- 
vations from 0.7 to 9 keV by I Young et"a l. (2002) show 
the hard X-ray flux to be peaked at a location consistent 
with that of the radio core and unresolved (less than 0."4 
in size, determined by comparison with a model PSF). 
Higher energy INTEGRAL observations show em ission 
between 20 and 100 keV (jBeckmann et al.l 120061 ). al- 
though there may be some contamination from intra- 
cluster gas. This hard X-ray emission is likely due to 
accretion disk emission Comptonized in the AGN corona, 
although the UV/optical emission is obscured. 

In addition to the AGN activity, HST imaging has also 
revealed star formation in the central region of Cy gnus 
A, which began < 1 Gyr ago (j Jackson et al.l [19981 ) . It 
is located in a 4 kpc ring around the nucleus, oriented 
orthogonal to the radio axis. 

Based on adaptive optics observations showing a sec- 
ondar y point source near the nucleus, Canalizo et al. 
(2003) suggest that Cygnus A may be in the late stages of 
a merger event. This merger event and related accretion 
may be related to the current epoch of nuclear activity. 

A near- infrared spect rum presented by 
iBellamv fe Tadhunterl ([2004D showed complicated 
emission line properties suggesti ng an infalling molec- 
ular cloud, consistent with the iCanalizo et al.l ()2003[ ) 
picture of a minor merger. The H2 lines were seen in 
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Table 1 

Mid-infrared Emission Line Properties 



V (THz) 

Figure 1. Spitzer IRS spectrum for Cygnus A. Extracted in a 
20" circular aperture from spectral mapping data. 

several components, both redshifted and blueshifted 
relative to the systemic velocity, interpreted as emission 
from a rotating torus. The observed near-infrared line 
ratios are consistent with excitation by X-rays (from 
the AGN), while likely ruling out shocks as a possible 
excitation method. 

For a more detailed summary of Cygnus A, including 
properties of t he larger host galaxy and environment, see 
the review bv lCarilli fc Barthell dl99l . 

While the bulk of the infrared emission is likely re- 
lated to the presence of an AGN, the star formation can 
contribute to the infrared emission as well. Addition- 
ally, the bolometric luminosity of the AGN is uncertain. 
Estimates have been made using template spectral en- 
ergy distributions (SEDs) and X-ray observations. How- 
ever, a significant portion of the bolometric luminosity 
comes out in the infrared, via dust reprocessing of the 
UV/optical continuum. An accurate determination of 
the bolometric luminosity then requires an understand- 
ing of the contributions to the observed infrared lumi- 
nosity. 

In this paper, we present modeling of the SED of 
Cygnus A. To accomplish this we combined new infrared 
measurements obtained by t he Spitzer Space Te lescope's 
Infrared Spectrograph (IRS; We rner et al.ll2004l ) with ex- 
isting Spitz er observations wi th the Imaging Array Cam- 
era (IRAC; iFazio et all 120041 ) and th e Multiband Imag - 
ing Photometer for SIRTF (MIPS; iRieke et all 120041 ), 
and measurements of the radio core from the literature 
to construct a radio through infrared (~ 2 — 10 5 GHz; 
4 — 10 5 fim) SED of the inner regions of Cygnus A. The 
resulting SED was subsequently modeled using compo- 
nents intended to replicate the physical processes likely 
to produce the observed emission. Using the results of 
the fitting we have been able to decompose the infrared 
emission and determine the bolometric luminosity of the 
AGN in Cygnus A. This in turn also provides an estimate 
of the star formation rate. 

The paper is organized as follows: in Section[2j we dis- 
cuss the new Spitzer IRS observations, new reductions 
of Spitzer IRAC and MIPS data, as well as the data 
from the literature. In Section [31 we describe the com- 
ponents used to model the SED, for which the results 
are described in Section |4j We conclude with a general 
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0.09 
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0.11 
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0.128 ±0.038 


0.31 




Note. — Only detections > 3tr are 


listed. Errors 



quoted from PAHFIT output. 

discussion in Section [5] and a summary in Section [6l 
2. DATA 

New observations were combined with data compiled 
from multiple archives and new reductions in order to 
obtain a SED covering almost five orders of magnitude 
in frequency. Our analysis is focused on the properties 
of the continuum emission in this system. Below we dis- 
cuss the reduction and analysis of the new Spitzer IRS 
observations and the re-reduction of archival IRAC and 
MIPS observations 

2.1. Spitzer IRS Observations 

Cygnus A was observed (PI: Baum) in "mapping 
mode" using the low resolution mode of IRS on board 
the Spitzer Space Telescope. Both the short- and long- 
wavelength slits were stepped across the source in incre- 
ments of half the slit width. After pipeline calibration at 
the Spitzer Science Center, the observations were com- 
bined into a spectral data cub e using the Cube Bui lder for 
IRS Spectra Maps (CUBISM: [Smith et al.|[2007al) . From 
this data cube, a spectrum was extracted in a 20" aper- 
ture (Figured]). 

2.1.1. Removal of Mid-infrared Emission Lines 
The IRS spectrum was fit using PAHFIT (|Smith et all 



l2007bl ) to measure and remove contributions from nar- 
row emission lines. Integrated line fluxes and widths are 
provided in Table [TJ While a study of the emission line 
properties is beyond the scope of this paper, we note the 
detection of multiple high ionization lines such as [O IV] , 
[Ne V] and [Ne VI I. These are all cons i stent with the pres- 
ence of an AGN ()Genzel et al1 H998: Ar mus et alJl2007L 
and references therein). Of particular note are [Ne V] 
and [Ne VI] which require the presence of ionizing pho- 
tons of at least 97.1 and 125.8 eV respectively. All the 
measured emission lines in Table [1] were subtracted from 
the IRS spectrum before fitting the SED. 

2.1.2. Dust Features 

In addition to fitting the emission lines, we have also 
fit the mid-infrared dust features using PAHFIT. Table 
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Table 4 

Flux densities from the literature 



Errors quoted from PAHFIT output. 



Table 3 

Additional Infrared Flux Densities from Spitzer 



Instrument / Channel 
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Aperture (") 


IRAC-2 


4.5 


0.010 


0.003 


12.2 


IRAC-4 


8 


0.054 


0.013 


12.2 


MIPS-70 


70 


2.20 


0.11 


30 


MIPS- 160 


160 


0.668 


0.033 


48 



[2] shows the integrated flux and profile FWHM based on 
the PAHFIT output. In contrast with the emission lines, 
the dust features were not removed as they are expressed 
i n the starburst mod els . 

ISpoon et al.l (|2007D developed a diagnostic diagram us- 
ing the 6.2 (an polycyclic aromatic hydrocarbon (PAH) 
and the strength of the 9.7 /zm silicate feature (S s y), 
which is seen in absorption in Cygnus A. These two spec- 
tral features can be used in tandem to classify the relative 
dominance of PAH emission, continuum emission, and 
silicate absorption. We follow their method of spline fit- 
ting to measure the depth of the silicate feature, finding 
S s iz ~ —0.8, where S s u is the negative of the apparent op- 
tical depth. The EQW of the 6.2 /im PAH is 0.0552 /im, 
placing Cygnus A on the border of region 1A and 2A in 
their diagnostic diagram, corresponding to objects dom- 
inated by continuum emission in the mid-infrared. This 
is consistent with the presence of a strong AGN. 

2.2. Archival Spitzer IRAC+MIPS Data 

The nucleus and ho tspots were imaged u sing IRAC at 
4.5 and 8.0 /mi (see [St awarz eTal1l2007l for a presen- 
tation of the data and analysis of hotspot properties). 
As fluxes for the core were not presented, the archival 
data were retrieved, re-reduced and calibrated according 
to the IRAC instrument manual. The images showed 
strong emission at the location of both radio hotspots 
and the radio core in both channels. The emission was 
unresolved in the core component, the measured flux den- 
sities for this component are given in Table[3] Extraction 
a pertures o f 12". 2 were used for both channels. 

iShi et all (|2005l ) presented observation of Cygnus A us- 
ing the MIPS instrument on Spitzer at 24, 70, and 160 
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/im. Their 24 /im flux is consistent with our IRS obser- 
vations. The slope between their 70 and 160 /im points is 
steeper than that of the Rayleigh- Jeans tail. The 70 and 
160 /im data were re-reduced from the BCD products in 
the Spitzer Science Center archive. Our measured flux 
densities are given in table [3] The core of Cygnus A was 
unresolved in both the 70 and 160 /im bands. Apertures 
of 30" and 48" were used at 70 and 160 /im, respectively. 
Our 160 /im measurement has larger error bars, but is 
broadly consistent with their value. 

2.3. Published Data 

Cygnus A's strong synchrotron emission at radio fre- 
quencies suggests that synchrotron may contribute to 
the infrared as well. To anchor the synchrotron spec- 
trum we supplemented the Spitzer observations with 
radio measurements f rom the lit er ature . Core fluxes 
were o btained from lEales et al.l (11989D: ISalter et al. 
(1989); lAlexander et al.l (|1984D ; iWright fc Birkinshaw 
(1984 ). Submillimeter n uclear fluxes were also taken 
from iRobson et all (|1998t ). Based on these archival data 
the unresolved radio cor^B is flat spectrum (a = 0.18, 
F v (x v^ a ) 1 until ~ 1 THz, where thermal emission from 
dust begins to dominate the SED. The compiled radio 
through infrared SED is shown in Figure [2] and the flux 
densities taken from the literature are provided in Table 
|U The resolution of the observations is also provided. 

There is some uncertainty associated with the highest 
frequency sub-mm observations. The 450 /im observation 
suggests that the synchrotron break may be occurring in 
this spectral regime. It is unclear if this is a genuine 
break or if the measurements are affected by variability. 
Additional observations at these frequencies may be able 
to clarify this issue. 

A consideration of the apertures is important when 
assembling data across several orders of magnitude in 
frequency. The resolution of the data used to assembled 
the SED in Figure [2] varies by over a factor of 10, but in 
all cases the core component is unresolved. At lower fre- 
quencies the emission is due solely to synchrotron emis- 
sion from the flat spectrum radio core. There is no evi- 
dence to suggest that the larger scale steep spectrum jet 
will contribute emission in the infrared. Accordingly we 
use the unresolved VLA-scale core to anchor the syn- 

9 Here "radio core" refers to the core seen by the VIA with on 
kpc-scale resolution. This encompasses flux from the VLBI scale 
core and jet. 
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Figure 2. Radio through mid-infrared SED with new Spitzer ob- 
servations and other data from the literature. See the text for 
details and references. 

chrotron emission in the nucleus, so the difference in 
apertures should not affect the construction of this SED 
for the nucleus. 

3. MODELING 

We aim to reproduce the major features in the SED: 
the powerlaw emission at radio wavelengths, the strong 
thermal emission at infrared wavelengths, and the over- 
all character of the silicate absorption. Extrapolating the 
powerlaw from the radio to the infrared requires a modi- 
fication of the powerlaw spectrum to avoid exceeding the 
observed mid-infrared flux. The infrared emission is ther- 
mal in nature, coming from dust at a variety of tempera- 
tures ranging from the cold ISM (T ~ 20 K) through hot 
dust near the AGN, up to the sublimation temperature 
(T ~ 1500 K). The power source for this dust heating is 
a combination of star formation and AGN activity, with 
an uncertain balance between the two. 

As noted in the introduction previous studies of 
Cygnus A have found evidence for simultaneous ongo- 
ing AGN activity and star formation, both of which can 
contribute to the infrared emission. We model the con- 
tinuum emission from ~ 2 — 10 5 GHz (3 x 10 6 — 5 /im; see 
Figure [2]) using components to represent the AGN torus 
model, a starburst, and synchrotron emission. The selec- 
tion of models has 15 free parameters, and the particular 
choices are discussed in the following sub-sections. 

3.1. AGN Torus Model 

According to the unified scheme for radio loud AGN, 
the SMBH and accretion disk can be hidden from view 
along some lines of sight by an obscuring torus. Along 
these obscured lines of sight the UV/optical radiation is 
absorbed and re-radiated in the mid- and far-infrared. 
As noted above, Cygnus A shows evidence for a hidden 
BLR through observations of Ha in polarized light. This 
suggests that Cygnus A harbors a "hidden" accretion 
d isk and a BLR which i s obscured along our line of sight. 

iNenkova et al.l ()2008l ) have constructed a model for an 
obscuring AGN torus, where the obscuration is due to the 
presence of multiple clouds along the line of sight to the 
AGN. We select this set of models because clumpy mod- 
els seem to provide better fits to Sil features than smooth 



dust distributions (e.g. B aum et al.l |2010T) . In order to 
reproduce the observed ratio of Type I /LI AGNs, the 
obscuring clouds collectively populate a rough toroidal 
structure with some opening angle. 

The CLUMPY torus model is specified by multiple ge- 
ometrical parameters. The outer radius of the torus (i? ) 
is Y times the inner radius (Rd), where the inner radius 
is determined from a dust sublimation temperature of 
T = 1500 K (R d = 0.4 x L ^ pc. L 45 is the bolometric 
luminosity of the AGN in units of 10 45 erg s _1 ). The 
geometry of the model is shown in Figure [3l The clumps 
have a Gaussian angular distribution, with a parame- 
terizing the width of the angular distribution from the 
mid-plane. The radial distribution is a power law with 
index q: r~ q . The inclination of the torus symmetry axis 
to the line of sight is i and the average number of clouds 
along a given line of sight is N. 
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Figure 3. Geometry of the CLUMPY torus model. Figure from 
Nenkova et al. (2008), used with permission. See the text for an 
explanation of the labels. 

In addition to the geometrical parameters, the mod- 
els also vary the optical depth of each clump (ry), us- 
ing the lOssenkopf et al.l (|1992 ) dust composition and 
IMathis et al.1 (|1977T) grain size distribution. The overall 
scaling of the model flux is Fagn, the bolometric flux of 
the AGN's accretion disk component (treating the syn- 
chrotron emission separately). 

The clumpy nature of the obscuration implies that 
there is a finite probability of a direct line of sight to 
the BLR, even at high inclinations. However, as Cygnus 
A shows no evidence for directly observed broad lines, we 
have only fit models where the central regions are fully 
obscured along our line of sight. 



3.2. Starburst 



Star 



the 
assume 



formation i s rep resented using 
ISiebenmorgen fc Kriigell ([2007) models which 
spherical symmetry and an ISM with dust properties 
characteristic of the Milky Way. Emission is broken 
down into two components: an old stellar population 
uniformly distributed through the volume and hot, 
luminous O and B stars embedded in dusty hot spots. 
The density of OB stars is centrally peaked although the 
model output is the emission integrated over the entire 
starburst. 
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The free parameters for the model are: starburst radius 
r, total luminosity Lsb, ratio of luminosity in O and B 
stars to total luminosity Job , visual extinction from the 
center to the edge of the nucleus Ay, and dust density 
in hotspots around O and B stars n. 

Owing to the spherical symmetry of this starburst 
model and the known ring morphology of the star for- 
mation in Cygnus A, parameters such as Ay and the ra- 
dius r do not have straight forward interpretations here. 
The assumption of optically thick star formation does 
not have strong evidence (either for or against), given 
the absence of high resolution far infrared observations. 

3.3. Synchrotron 

The strong radio emission in Cygnus A justifies the 
final model component. VLA core fluxes are consis- 
tent with powerlaw emission to the submillimeter where 
thermal emission begins to dominate. Extrapolating the 
powerlaw to higher frequencies suggests that the syn- 
chrotron spectrum must either be modified or subjected 
to attenuation. The AGN is known to sit behind a 
dust lane with signifi cant extinction (Ay = 50 ± 30; 
iDjorgovski eTaO l991). Though some attenuation of the 
synchrotron flux is expected at higher frequencies, this 
alone is not sufficient to explain why the synchrotron 
powerlaw does not continue through the mid-infrared. 
With extinction alone, the flux densities at 30 and 10 n 
m would be similar for the observed spectral index of 
a = 0.18. For Ay — 50, the expected flux from the ex- 
trapolated synchrotron emission would by itself exceed 
the measured infrared flux at 10/zm. We conclude that 
there must be a break in the population of emitting elec- 
trons. 

One possible explanation is a simple cutoff in the pop- 
ulation of relativistic electrons at some energy (Case I). 
This would manifest itself as a cutoff in the synchrotron 
spectrum: 

f,oci/ _ai r^e^M (1) 

where v c is the frequency corresponding to the cutoff in 
the energy distribution of the particles, a.\ is the spectral 
index in the optically thin regime, and r r (y) is the dust 
screen between the sy nchrotron emitting re gion and the 
observer. r r {y) follows iDraine fc Led (|1984j) . 

In Cygnus A the spectral index ai is measured from 
VLA core radio fluxes, and the amplitude of the pow- 
erlaw fixed from the same observations. The only free 
parameters are the cutoff frequency v c and the optical 
depth r r . These parameters art partly degenerate in that 
in trial runs we experienced a situation where t t and v c 
would both increase, effectively offsetting each other. To 
combat this we limited r T such that Ay does not exceed 
500 towards the radio source. 

An alternate model for the synchrotron emission at 
higher frequencies is a broken powerlaw behind a dust 
screen (Case II). Aging of the population of relativistic 
electrons results in a broken powerlaw whose spectral 
index in creases for frequ encies higher than a break fre- 
quency (|Kardashevlll962D . The functional form adopted 
is: 

\V 012 if V > V break 



For a simple aging of the electron population, the post- 
break spectral index in Cygnus A would be cti — 1.24. 
As with Case I, we fit a dust screen in front of the syn- 
chrotron emission (r r ). 

3.4. Stellar Contribution to the Mid-infrared 

Starlight can potentially contribute to the mid-infrared 
flux, especially in a larg e aperture. Using the flux of the 
stellar component from Jackson et al. ( 1998) and their 
brightness profile, we determined the contribution due to 
starlight in a 20" aperture. The relative fluxe s were con- 
sistent with an elliptical galaxy template from lSilva et al.l 
(1998). Expected flux densities at 2.2, 5, and 10 /Ltm were 
computed using the same template spectrum. After sub- 
tracting a nuclear point source, the K-band flux density 
in starlight is consistent with the flux in a 20" aperture 
as measured us ing Two Micron All Sk y Survey (2MASS) 
data products (Skruts kle et al.1 12006). At 5 and 10 pm, 
the starlight contributes roughly 14% and 2% respec- 
tively, of the flux measured by IRS. Therefore, we do not 
include any contribution to the mid-infrared flux from an 
old stellar population in our modeling. 

3.5. Dust in the NLR 

Warm dust in the NLR directly illuminated by 
the AGN can also contribut e to the m id - an d 
far-infrared emission (e.g ., IGroves et al.l l2006h . 
I Ramo s Almei da et al.l (|2009T ) find that a significant 
contribution to the mid-infrared flux can come from 
sources other than the AG N torus, particularly addi- 
tional hot dust. In Cygnus A lRadomski et al.l (|2002ft find 
an extended component to the mid-infrared emission 
which is consistent with T ~ 150 K dust. Some of this 
extended emission is co-spatial with sites of possible star 
formation. Dust of this temperature is well reproduced 
with our choice of starburst models (with the implicit 
assumption that this dust is heated by star formation). 

Higher temperature dust is also present in the inner 
regions of Cygnus A. However t he resolved mid-infrared 
images of the nuclear regions bv lRadomski et al.l ((2002) 
are unable to distinguish between dust in the "torus" and 
dust in the NLR heated by the AGN. Without strong ob- 
servational motivation for an additional dust component 
we opt not to include one. As will be shown in the next 
section we are able to reproduce the observed emission 
without this additional component. 

3.6. Prior Constraints on Model Parameters 

Prior to fitting the models, the available parameter 
space was constrained using results from previous studies 
of Cygnus A. The opening angle and covering fraction 
of the CLUMPY torus are determined from the a and 
N pa rameters (see Equations (3) and (4) in iMor et all 
2009). We define the half opening angle of the torus as 
the angle at which the escape probability of a photon 
drops below e : 

Ohaif = 90-a^hJf (3) 

iTadhunter et al.l ((19990 measured the opening angle of 
the ionization cone in Cygnus A using near infrared HST 
data, finding Ohaif — (58 ± 4)°. We use this value as the 
torus opening angle and limit the parameter range of a 
and N according to Equation [3] 



6 



100000 10000 



Wavelength (^im) 
1000 100 




v (THz) 



10 



Wavelength (|am) 
100000 10000 1000 100 



0.1 - 



0.01 



0.001 - 



0.0001 



0.001 




v (THz) 



1000 



Figure 4. LEFT: Case I fit to the Cygnus A SED. The shaded region shows all acceptable model fits within the 95.4% confidence interval. 
The lines are the components for the best-fit (lowest X 2 ). RIGHT: Same as left, but for Case II fit to the Cygnus A SED. Insets are a 
zoom of the region around the 9.7 /im Sil absorption. 
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Figure 5. Histogram of bolometric AGN luminosity from best fit torus parameters (left) and bolometric starburst luminosity (right) for 
Cases I (blue, solid line) and II (red, dashed line). 



The inclination range of the torus i was limited by 
very long baseline interferometry (VLBI) observations 
and modeling of the inner pc-scale jet (50° < i < 85°, 
iSorathia et al.lfl996l ). Additionally, the available VLA 
core radio fluxes were used to fix the synchrotron spec- 
trum at radio frequencies (ai, and the amplitude). 

4. RESULTS 

To fit the observed SED, the flux from each model was 
summed in each wavelength bin and compared to the ob- 
served flux within a 20" aperture. Optimization was per- 
formed using Levernberg-Marquardt least-squares mini- 
mization. 

Some degeneracy between parameters was seen in the 
model results. Comparison model runs with fewer prior 
constraints on the parameters resulted in similar % 2 val- 
ues to those quoted below for regions of parameter space 
which are unlikely to be physically reasonable matches 
to Cygnus A (e.g., CLUMPY torus fits with covering 
fractions of unity and O^aif = 0°). Prior constraints 
from other observations are thus important in eliminat- 
ing regions of parameter space which may be statistically 
reasonable but physically unrealistic. 

After the models were fit, confidence intervals were de- 
termined separately for Case I and Case II through boot- 
strapping with replacement (|Efronl ]l981l ). These were 



used to determine the "acceptable" range of parameter 
values. 

The results of the fits are shown in Figured] The range 
of parameter values for fits within at 95.4% confidence 
interval are shown in Figures [5]|9] Section 14.11 discusses 
the Case I fits (exponential cutoff in the relativistic elec- 
tron population) and Section 14.21 discusses the Case II 
fits (broken power law for the synchrotron emission). 

4.1. Case I: Synchrotron Exponential Cutoff 

The best fit combination of parameters for Case I has 
X 2 /DOF = 1.35. The range of acceptable matches for 
the 1641 model combinations is shown in Figured] (left). 
Histograms of the best fit parameters are given in Figures 
[5] (left) and[Hl where / is the fraction of models within the 
95.4% confidence interval in a given bin. Parameters for 
the best fit model are marked with a small blue horizontal 
bar. 

4.1.1. AGN '/Torus Properties and Contribution 

The fits favor a bolometric accretion disk luminosity 
of log(L AGN /L Q ) ~ 11.8 - 12.0 (median of 11.82 with 
an interquartile spread of 0.09). The fits clearly favor 
an extended torus (Y — 200) with low opacity clouds 
(Ay ~ 10 — 30). For the median bolometric luminosity, 
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Figure 6. Histograms of torus parameters for acceptable fits within the 95.4% confidence interval for Cases I (blue, solid line) and II 
(red, dashed line). The small horizontal bars denote the parameter value for the best-fit model. From left to right, top to bottom: viewing 
angle of the torus (90 meaning the torus is viewed edge-on), N - average number of clouds along an equatorial line of sight, Ay - for an 
individual cloud, a - angular width of the distribution of torus clouds, Y - ratio of inner to outer radii, q - power law index of the radial 
distribution of torus clouds, covering fraction of the torus (equivalent to the escape probability of optical/U V photons, assuming optically 
thick clouds), integrated Ay along an equatorial line of sight, integrated Ay along our line of sight (calculated using N, a, the cloud Av, 
and the inclination). Parameters in the first two rows are discrete and the histograms represent the relative number of models with those 
specific values. Parameters in the bottom row are derived from the parameters in the top two rows as well as the luminosity of the AGN 
component (shown in Figure [5}, and are continuous. 

the inner radius of the torus is Rd — 0.6 pc, giving a 
corresponding outer torus edge of ~ 125 pc (for Y= 200). 

The range of values for i was limited based on previous 
work in the radio regime. The fits prefer an inclination 
for the torus on the high end of the range, i ~ 80°. 
N has a bimodal distribution, with 25% of fits having 
N < 6 clouds along an equatorial line of sight and ~ 
60% having N > 20 clouds. Most of the fits within the 
95.4% confidence interval have a flat radial distribution 
of clouds (q = 0). 

The opacity of individual clouds is anti-correlated with 
the radial extent of the torus; small torus sizes favor high 
Ay values. High values of N (i.e., many clouds along 
a given line of sight) are favored for small torus sizes, 
however for large torus sizes, both small and large N 
values are acceptable. 

We calculate the total extinction through the torus 
along both an equatorial line of sight (Figure [6] middle, 
bottom row), and along our line of sight to the torus 
(Figure fright, bottom row). The integrated equatorial 
Ay spans a range between 100 and 250. When view- 
ing angle is taken into consideration, the range narrows, 
with most fits showing the line of sight Ay between 80 



and 120. However the Ay is poorly constrained beyond 
having a well-defined lower bound. 

The tor us covering fraction was computed using the 
method of lMor et alT(|2009l their Equations (3) and (4)). 
By design, the derived covering fractions for the torus are 
between 50% and 70%. 

4.1.2. Starburst 

The median value of the starburst luminosity for Case I 
is \og{LsB / Lq) ~ 10.8, with a tail up to ~ 11.6 contain- 
ing approximately 40% of the fits (Figure [5]). This covers 
a range of star formation rates from 10 to 70 yr " 1 



as dete rmined by the Lir calibration from 
(I1998bl) (for a review of SFR estimates, see 
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Kennicutt 
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Roughly half the fits within this confidence interval 
have 40% of the starburst luminosity in the form of 
OB stars. The distribution of size parameters is rela- 
tively flat, although even the largest sizes would be un- 
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Figure 7. Histograms of starburst parameters for acceptable fits within the 95.4% confidence interval for Cases I (blue, solid line) and 
II (red, dashed line). The small horizontal bars denote the parameter value for the best-fit model. From left to right, top to bottom: the 
fraction of luminosity in O and B stars, integrated extinction through the starburst, size of the starburst, dust density in hotspots around 
O and B stars, fo s and the size are discrete and the histograms represent the relative frequency of models with those specific values. The 
Ay and n parameters are more finely sampled than the histogram represents; the plots show the relative fraction of models with parameters 
in the noted ranges. 



resolved by our Spitzer observations. The distribution of 
dust density n peaks around 10 3 cm' 3 . The extinction 
through this starburst is relatively unconstrained. 

The strengths of dust features in starburst models 
falling within the 95.4% confidence interval were also 
measured using PAHFIT to compare with the intensi- 
ties of observed dust feature. Figure |8] shows histograms 
of the predicted intensities from our models divided by 
the measured intensity from the IRS spectrum. Only 
dust features which are detected in the IRS spectrum 
have been plotted. In general the agreement is good 
with most models matching the measured line intensi- 
ties. However for some dust features a significant num- 
ber of the models within the confidence interval predict 
emission in excess of what is observed. While a detailed 
comparison of the relative strengths of dust features is 
beyond the scope of the paper, we suggest the models 
are able to generally reproduce the dust features seen in 
the spectrum. The dust features which show the greatest 
discrepancy between observed fluxes and those predicted 
from the modeling (6.2, 7.8, and 11.3 /xm) comprise 3 of 
the 5 lowest signal-to-noise dust feature fits in the spec- 
trum. The discrepancy may then be due to the difficulty 
of fitting low equivalent width dust features. 

PAHFIT also attempted to fit eight other dust fea- 
tures in the IRS spectrum. The upper limits for three 
(5.7, 14.0, and 15.9 /im) are consistent with expectations 
from the starburst models. One dust feature (16.4 /im) 
has a measured limit below that expected from the best 
fitting starburst model, and so is discrepant. The other 
four features (7.6, 10,7, 11.2, and 12.7 fim) also have 
upper limits from PAHFIT which are lower than the ex- 
pected value from the starburst models. However, these 
are coincident with or near other spectral features (e.g., 
narrow emission lines or the Sil absorption). Thus an 
accurate measurement of these dust features in the IRS 



spectrum would rely more heavily on fitting the wings, 
which could prove difficult in a continuum dominated 
source such as Cygnus A. 

4.1.3. Synchrotron Radiation 

The synchrotron emission amplitude was fixed using 
the non-thermal emission from the radio core, assuming 
it to be a point source at all frequencies observed. The 
synchrotron model, therefore, has only two parameters: 
cutoff frequency (y c ) and extinction due to dust (r r ). 

The distribution of cutoff frequencies is somewhat 
broad, covering the range of v ps 10 — 60 THz (5 — 30 /im) 
(Figure[S]) . The fits show a range of acceptable extinction 
for the dust screen, peaking in the range of Ay w 60 — 80, 
slightly lower than the predicted line of sight Ay from the 
torus models. The integrated luminosity of the "core" 
synchrotron flux is log(L sync / Lq) ~ 11.2. 

This break in the synchrotron spectrum is consistent 
other powerful FR II radio sources, where an extrap- 
olation of the radio synchrotron emissi on significantly 
exceeds the observed optical flux (e.g., Schwartz et al.l 

2009]). In these 



2000; Sa mbruna et~alT 2004; M ehta et all 



cases, the similarity of the radio and X-ray spectral in- 
dices suggests that the X-ray may be produced by inverse 
Compton (IC) scattering. 

With the modeled synchrotron spectrum, a magnetic 
field strength, and the assumption that each electron 
emits at a single frequency, the energy distribution of 
relativistic electrons can be determined: 



7(1/) 



47rm e ci^ 
3eB 



(5) 



where all constants and values are in the cgs system of 
units. 
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The magnetic field strength has been estimated from 
VLB1 obse rvations to be rough l y 17 mG (100 mG) in th e 
jet (core) (jRoland et al.l 119881 iKellermann et aLlll981[ ). 
Using v ss 30 THz and Equation [5] we find 7 - 2 x 10 4 
(8 x 10 3 ) for the jet (core). 

In addition to radiating energy via synchrotron emis- 
sion, relativistic electrons can also lose energy via IC 
scattering. The ratio of the energy lost through syn- 
chrotron to the energy lost to IC is simply given by 
the ratio of the magnetic and radiation energy densities, 
with the cosmic microwave background (CMB) setting 
a minimum value for the radiation energy density. Us- 
ing the magnetic field strength assumed above (17 mG), 
and Tcmb = 2.73 K, synchrotron losses are dominant 
by a factor > 10 8 . A larger magnetic field would fur- 
ther enhance this ratio, so IC losses from scattering of 
CMB photons have a negligible effect on the population 
of relativistic particles. 

4.2. Case II: Synchrotron Broken Power Law 

The range of best fit models for Case II is shown in 
Figure H (right). The best fit model had x 2 / D0F = 
1.03. Figure 0] (right) shows the range of the 1580 model 
fits within the 95.4% confidence interval. Histograms of 
the best fit parameters are given in Figures [5] (left) and 
[6l A red horizontal bar marks the parameter for the 
parameters with the lowest % 2 value. 

4.2.1. AGN/Torus 

In Case II fits the bolometric accretion disk luminos- 
ity is roughly \og(L,AGN / Lq) ~ 11.8 (median 11.88 with 
an interquartile range of 0.08). The distribution of lu- 
minosities is more strongly peaked than for Case I, with 
over 80% of fits occupying the log(L/L ) = 11.8 - 12.0 
bin. A large torus (Y = 200) is exclusively preferred, 
corresponding to R ou t ~ 135 pc (using the median value 
for the luminosity). Again, by design, the torus covering 
fraction is between 0.5 and 0.7. An inclination of i — 80° 
is exclusively preferred. 

The distribution of the average number of clumps along 
a line of sight N is again bimodal with either large (< 
14) or small (< 8) numbers of clouds preferred. The 
number of clouds is anti-correlated with the extinction 
through individual clouds. The spread in equatorial Ay 
is somewhat large (120 — 260), but the line of sight Ay 
is more constrained (100 — 160 for 80% of fits). 

4.2.2. Starburst 

The starburst component in Case II fits show qual- 
itatively similar behavior to Case I. The typical lumi- 
nosity is comparable, but shows a smaller tail up to 
log(Ls_B/i©) = 11-6 (~ 30% of fits). Other parame- 
ters are similar to those in Case I fits, with the exception 
of the size which appears to have a very slight preference 
for a larger size, suggesting an overall cooler dust tem- 
perature for a given luminosity. This may indicate a de- 
generacy between the starburst and synchrotron models. 
The critical frequency for the synchrotron spectrum (see 
the next section) is at lower frequency for Case II, result- 
ing in a smaller contribution to the far-infrared flux. The 
starburst model compensates with a larger overall size to 
provide additional cool dust emission (for a given star- 
burst luminosity). The dust feature strengths expressed 
in the models generally compare favorably with Case I. 



4.2.3. Synchrotron Properties and Contribution 

An alternate mechanism for limiting the influence of 
synchrotron emission at shorter wavelengths is for the 
spectrum to break at some frequency (see Section l3~3f . 
Case II fits used a fixed pre-break spectral index of 
a.\ = 0.18 and a post-break spectral index of cx-i = 1.24, 
consistent with aging of the relativisti c population (with - 
out injection of additional particles; iKardasrievl [1962) . 
As the flux density decreases in a slower fashion when 
compared to an exponential cutoff, the powerlaw must 
break at lower frequencies to ensure that the observed 
mid-infrared flux is not exceeded. 

The break frequency is roughly 5 THz (60 fim). 
The unobscured synchrotron luminosity is found to be 
\og(L sync / Lq) ~ 11.1. Following the same arguments as 
Case I, the electrons emitting at the break frequency have 
7~ 7 x 10 3 (3 x 10 3 ) for jet (core) magnetic field values. 
Again, this is consistent with the observed properties of 
jets in other FR II radio sources. 

The extinction of the dust screen in front of the ra- 
dio source is lower than Case I fits, with Ay < 60, still 
within ran ge of the extinction to t he central source es- 
timated bv iDiorgovski et all (|1991f ). but lower than the 
computed line of sight Ay from the torus. The larger dis- 
crepancy compared with Case I is due to the enhanced 
short wavelength emission of the broken power law at 
shorter wavelengths (when compared to the Case I expo- 
nential cutoff). The equatorial torus Ay increases to pro- 
vide an overall cooler temperature, and the modeled dust 
screen in front of the synchrotron component remains 
small (in order to contribute sufficient flux at shorter 
wavelengths). 

5. DISCUSSION OF GENERAL SED RESULTS 

5.1. Luminosity and Kinetic Power in Cygnus A 

Our best estimate of the bolometric AGN luminos- 
ity in Cygnus A (log(i/L ) ~ 12, including the syn- 
chrotron component a n d X-r ay emission) is above the 
iWhvsong fc Antonuccil (|2004f) estimate using Keck mid- 
infrared observations and a PG-quasar spectrum (3.9 x 
10 11 Lq), but is somewhat below the estimates of 
iTadhunter etafl (|200l who find L bol = 12 - 55 x 10 11 
Lq (although consistent with the low end of their range) . 
The bolometric AGN luminosity is the best constrained 
parameter, and is insensitive to the synchrotron model 
adopted. 

The kinetic power in the expansion of the radio lobes 
in Cygnus A can be infe rred from X-ray obs ervations of 
the cluster environment (| Wilson et al.l[2006h . From this 
analysis of Chandra observations Cygnus A has a kinetic 
power of L kin = 1.2xl0 46 ergs" 1 (log(L fein /L ) = 12.5), 
a factor of three larger than the bolometric AGN lumi- 
nosity inferred from the modeling. Given the uncertain- 
ties in the determinations of both the kinetic and bolo- 
metric power, it is unclear if the kinetic power dominates 
significantly over the power emitted as radiation. 

5.2. Implications of the Modeling 

The probable torus sizes from the SED modeling give 
outer radii of R a = 130 pc (« 0".2 at the distance of 
Cygnus A) . Generally, the torus parameters provide rea- 
sonable estimates of the properties of the obscuring struc- 
ture. This predicted outer radius is significantly larger 
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than what is generally assumed to be reasonable for the 
obscuring torus. A torus with the angular diameter of 
the w 0."2 is well within range of multiple mid and near- 
infrared instruments. At large radii, the distinction be- 
tween torus clouds and narrow line regio n clouds may be 
somew hat arbitrary. The observations bv lCanalizo et aT] 
(2003) show extended emission on scales including and 
larger than our fit torus sizes. The degree of contamina- 
tion by emission lines is unclear, and it would be inter- 
esting to attempt to obtain a line-free continuum image 
to ascertain the true size and shape of the continuum 
emitting region. 

The torus size parameters are consistent with those 
found from modeling of the 9.7 fim Sil feature by 
llmanishi fc Uenol ()2000D . Their radiative transfer mod- 
eling suggested some torus properties similar to those 
presented here, namely a small inner radius (< 10 pc) 
and an inner-to-outer radius of 80 — 500. In contrast 
they find a steeper radial dependence for the dust distri- 
bution q ~ 2 — 2.5, possibly due to the use of a smooth 
dust distribution. 

The preferred parameters for the obscuring torus; low 
q, high Ay, and large Y suggest the model is being driven 
toward a cooler spectrum. This could influence the star- 
burst component, leading to a lower inferred star forma- 
tion rate. In order to combat this it would be beneficial 
to place tighter constraints on the far-infrared SED to 
enable better modeling of the star formation in Cygnus 
A. 

A comparis on of our IRS d a ta wit h mid-infrared ob- 
servations by rRadomsk i et al.l ((2002) is consistent with 
the suggestion of our model that at 10 and 18 /j,m, the 
emission is dominated by the t orus+synchrotron com- 
ponent. iTadhunter et al I (|1999t ) measure a K-band nu- 
clear point source with F v — (49 ± 10) /zJy. This flux 
limit is broadly consistent with some model fits for the 
torus+synchrotron component. 

For the estimated bolometric luminosity 10 12 Lq, 
and a black hole mass of 2.5 x 10 9 M , the Edding- 
ton Ratio (L/L edd ) is - 1.3 X 10~ 2 . This is simi- 
lar to, but sligh t ly low er than previous estimates from 
ITadhunter et al.l (|2003f ). likely due to the fact that our 
model attributes some of the IR luminosity to the heating 
of dust from star formation. This starburst luminosity is 
relatively well constrained, with consistent values across 



the bulk of fits. The AGN (torus plus jet) contributes 
~ 90% of the infrared luminosity and star formation pro- 
duces the remaining 10%. 

5.3. Future Work and Observations 

The results of our modeling can be tested and improved 
with the help of future observations in various wavelength 
regimes. 

Around 1 THz (~ 300/xm) the synchrotron and star- 
burst model components are of similar flux density, with 
the synchrotron contribution decreasing and the thermal 
contribution from star formation beginning to dominate. 
Unfortunately, data at this location are limited in reso- 
lution, sensitivity, and wavelength coverage. Disentan- 
gling the emission from the starburst and emission from 
the AGN at this frequency will be possible with the Ata - 
cama Large Millimeter Array (ALMA; e.g., Carilli 2005), 
particularly with the availability of Band 10 and full sci- 
ence capabilities. 

The far-infrared suffers from poor sampling of the 
SED. In this regim e the Herschel Space Observatory 
(jPilbratt et alJl2010( ) provides an opportunity to improve 
on the understanding of the continuum emission. SPIRE 
and PACS cover wavelength ranges of interest: the re- 
gion of the possible jet-break and the long wavelength 
side of the infrared bump. Observations here can provide 
improve d constraints for the synchrotr on and starburst 
models. iGonzalez- Alfonso et al.l (|2010l ) demonstrate de- 
composition of far- infrared emission for Mrk 231, with 
the Herschel SPIRE observations providing important 
constraints. Similar observations of Cygnus A will pro- 
vide important data on this sparsely observed portion of 
the SED. PACS observations would contribute measure- 
ments of far-infrared fine structure lines which could be 
used in concert with the mid-infrared lines to provide an 
alternate method of determining the relative contribu- 
tion of star for mation and AGN act ivity to the infrared 
luminosity (e.g.. iFischer et al.|[20Tol , for Mrk 231). 

Although the synchrotron spectrum breaks in our fits, 
it may still dominate the flux between 5 and 10 /im. This 
result can be tested by future mid-infrared polarimetry 
or variability studies. Synchrotron emission from a com- 
pact source such as the radio core or jet knots should 
be subject to flaring. If the emission in this wavelength 
range is instead of a thermal origin (e.g., hot dust in the 
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host galaxy), the 10 /jm flux will remain relatively stable. 
6. SUMMARY 

Using a combination of a new mid-infrared spectrum 
from the Spitzer Space Telescope and radio data from 
the literature, the nuclear emission in Cygnus A has 
been modeled combination of powerlaw emission 

from a synchrotron jet, reprocessed AGN emission from 
a dusty torus, and emission from a dusty circumnuclear 
starburst. 

The data are well fit by a combination of these three 
models, and all three are necessary to reproduce the ob- 
served emission. Statistically acceptable fits were found 
for both an exponential cutoff in the population of rel- 
ativistic electrons (Case I) as well as emission from an 
aging electron population (Case II), however we are un- 
able to distinguish between the two cases. For Case I 
we find the cutoff frequency to be between 10 and 50 
THz (5 — 30 /Km) while for Case II fits the predicted 
break frequency is 5 THz (60 /im). Degeneracy between 
the starburst and synchrotron components makes a more 
precise determination of the break or cutoff frequency dif- 
ficult. Better observations on the long- wavelength side 
of the thermal bump will provide tighter constraints on 
the models. 

From this modeling, we find the following: 

• The bolomctric luminosity of the AGN in Cygnus 
A is - 10 12 L Q 

• The mid-infrared emission is consistent with emis- 
sion from a clumpy obscuring torus with an outer 
size of ~ 130 pc. 

• The far-infrared emission is consistent with being 
dominated by star formation which is occurring at 
a rate between 10 and 70 Mq yr _1 . 

In Cygnus A, the infrared emission is a combination 
of AGN and starburst heated dust, with the AGN con- 
tributing ~ 90% of the luminosity. 
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